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ABSTRACT 



We calculate the damping of quadrupole f- and low order g-modes (primary 
modes) by nonlinear coupling to other modes of the star. Primary modes 
destabilize high degree g-modes of half their frequency (daughter modes) 
by 3-mode coupling in radiative zones. For sunlike stars, the growth time 
= jy -1 pa 4i? 4^ days, where £0,42 is the initial energy of the primary mode 
in units of 10 42 erg, and the number of daughter modes N ~ 10 10 i?o,42- The 
growth rate is approximately equal to the angular frequency of the primary 

1 /2 

mode times its dimensionless radial amplitude, 5R/R* ps 0.002£ , 0j / 42 . Although 
the daughter modes are limited by their own nonlinearities, collectively they 
absorb most of the primary mode's energy after a time ~ 10?7 1 provided 
E > 10 40 erg. This is orders of magnitude smaller than usual radiative damping 
time. In fact nonlinear mode interaction may be the dominant damping process 
if Eq ^ 10 37 erg. These results have obvious application to tidally captured main 
sequence globular cluster stars of mass > 0.5Mq; the tidal energy is dissipated 
in the radiative core of the star in about a month, which is less than the initial 
orbital period. 

Nonlinear mode coupling is a less efficient damping process for fully 
convective stars, which lack g-modes. In convective stars most of the tidal 
energy is in the quadrupole f-modes which nonresonantly excite high order 
p-modes of degree 0, 2, and 4. The resultant short wavelength waves are more 
efficiently dissipated. The nonlinear damping time for f-modes is shown to be 
proportional to 1/E ; this damping time is about 30 days for E ps 10 45 erg 
expected in tidal captures. However, at such a large energy the system is very 
nonlinear: four-mode and higher-order couplings are as important as three-mode 
couplings. 
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1. Introduction 



Tidal interactions in highly eccentric close binary systems primarily excite low order 
quadrupole gravity modes (g-modes) and the fundamental mode (f-mode), because these 
modes have periods comparable to the periastron passage time and their eigenfunctions 
have the largest overlap with the tidal forcing field [[Kumar, Ao, fc Quataert (1995)|| . The 



dissipation time scale for these modes, for main sequence stars, are 10 4 years or greater. A 
long damping time typically leads to a larger mode amplitude for modes that are subjected 
to repeated tidal forcing, but perhaps more interestingly offers the possibility that in tidal 
interactions energy can be transferred from the modes to the orbit as well as from the orbit 
to the modes, causing a chaotic evolution of the orbit ||Mardling (1995) 



m some cases. 



In the case of tidal capture binaries, this process can lead to ionization of the system in 
subsequent encounters fKochanek (1992)| , |Mardling (1995)| . 



Given the large amplitudes of modes in close binary systems, one expects nonlinear 
interactions amongst modes carrying most of the energy and others, which are tidally 
not excited, to be an important process of draining energy from the primary modes. We 
investigate this process and assess its significance in close binary systems. 

The next section describes the basic equations for mode couplings and discusses the 
most significant nonlinear process for energy transfer from modes of a radiative star. In 
section §3 we discuss the nonlinear mode couplings for a convective star. Applications to 
tidal capture are discussed in §4, and our main results are summarized in §5. The details of 
the calculation of nonlinear coupling coefficients are presented in Appendix A. 



2. Nonlinear Interactions: Stars with Radiative Zones 



The leading-order nonlinear interactions of stellar modes are 3-mode couplings. These 
can be derived by expanding the Lagrangian density for adiabatic motions in powers 
of the fluid displacement. If the expansion is truncated at second order, the standard 
Euler-Lagrange procedure yields linear equations of motion, which one solves for the normal 
modes. When the third-order parts of the Lagrangian density are retained, the equations of 
motion for the mode amplitudes contain quadratic terms: 

A a + U 2 a A a = -U 2 a KafljApAy, (1) 
0,7 

where uj a & A a are mode frequency and amplitude. The 3-mode coupling coefficient, K a p-y, 
is completely symmetric in the three indices. It is obtained by substituting the normal 
mode expansion of fluid displacement into the cubic Lagrangian density, and integrating 
over the stellar interior (see Appendix A). Since our eigenfunctions are normalized to unit 
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energy, the energy in mode a is A 2 a . Note that equation ([I]) with its right side set to zero 
describes free normal modes of oscillation. We could represent tidal forcing by adding a 
term independent of the amplitudes, but in this paper we are concerned with the evolution 
of the modes after forcing has ceased. 

We assume that initially all the pulsation energy resides in a few low order modes. If 
the star has a radiative core, then these modes couple to high degree g-modes and transfer 
energy away by parametric resonance instability. This is discussed in the next sub-section. 
Even for completely convective stars nonlinear energy transfer, involving non-resonant 
3-mode interaction with high frequency p-modes, can be an important process and is 
discussed in 53. 



2.1. Parametric instability 



Parametric instability arises when a mode is nonlinearly coupled to another mode of 
approximately twice the frequency, or more generally when two low frequency modes of 
infinitesimal amplitude couple to a higher frequency mode of finite amplitude. To analyze 
this case we consider a single triplet interaction involving modes a, (3\ and /?2- Expressing 
mode amplitudes in the following complex notation 

A a = Aexp(iu a t) + c.c., Ap 1 = £>i exp^co^t) + c.c, Ap 2 = B 2 exp(iup 2 t) + c.c, (2) 

and substituting this into eq. (|1|) we obtain the following simplified equations 

Bi+jBx/2 = i3ujp 1 K a p l p 2 AB2exp(iAut), (3) 
^2 + 7-82/2 = i3up 2 Kai3 1 i3 2 ABlexp(iAut), (4) 

where Au = u a — ujp 1 — uip 2 , 7 is the linear dissipation rate of these modes, and the factor 
of 3 appears for combinatorial reasons. Rapidly oscillating terms as well as second time 
derivative of mode amplitudes were discarded in deriving these equations. At initial times 
when the amplitudes of modes (3\ and (3% are infinitesimal, we can take the amplitude of 
mode a to be fixed (A ), and solve the above two equations for the growth rate (77) of mode 
amplitudes, which is 

V = 2 yVAWahfhVPiVfh ~ ( AuJ ) 2 ~ 7] • ( 5 ) 

The coupling coefficient k Q/ 3i/3 2 is largest when the modes /3i and fa have nearly the same 
radial order n and spherical-harmonic degree £, whence Up 1 ~ Up 2 ~ ui a /2. We see that 
parametric instability sets in provided that 

(Auf < 9K 2 a p 1 p 2 A 2 Q ujp 1 ujp 2 - 7 2 , (6) 
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Ignoring the damping for now, we obtain a very simple criterion for parametric instability: 

|Au;| < 3u a A K al3l p 2 /2. (7) 

The product AqK^^ is a dimensionless quantity which is a measure of the nonlinearity 
of the system. The coupling coefficients are calculated using equation QA17) of Appendix 



A, which is a somewhat long expression. However, simple scaling argument shows that 
K ai3j3 ~ G~^ 2 rV 2 ''M^ 1 for a low order quadrupole mode coupled to two high order 
g- modes (where M* and i?* are stellar mass and radius), so that K a/3lf32 « 10~ 24 erg -1 / 2 
for a solar-mass main sequence star. This is in good agreement with detailed numerical 
calculations (Table 1). If the energy E = A$ in the f-mode is 10 42 erg then parametric 
resonance is possible provided that we can find a triplet with Au less than 0.15% of the 
frequency of the f-mode. We discuss below the possibility of finding such resonant triplets 
involving high degree g-modes. 

To put this energy in perspective, note that the dimensionless radial amplitude of the 
f-mode is 5R/R* ~ 2 x 10~ s Eq%, where E 042 = -Eo/(10 42 erg). Tidal capture is expected 
to result in motions 5R/R* ~ 0.1, E ~ 10 45 erg. For comparison, the most energetic 
five-minute oscillations of the Sun have about 10 28 erg, and the total energy in all solar 
modes is ~ 10 33 erg | |Goldreich, Murray, fc Kumar (1994)| |. It also follows from these 



scalings that Alu/lu ~ 2r]/ui ~ 5R/R 

The dispersion relation for g-modes is [e.g. phristensen-Dalsgaard fc Berthomieu 



199T) | 



2^£{£ + 1 

where 



Pni= r ^^ (2n + £-5) + ^- 



Vo 



+ 1 



2tt 2 



(9) 



u J^drN B y 

P n e (= 2ir/u n £ = l/vni) is the mode period, Nb is the Brunt- Vaisala frequency, r\ is the 
upper turning point of the mode, 5 ~ 5/6, and V\ and V 2 are constants of order 0.5 and 
6 respectively for 1 Mq main sequence star. Note that mode frequency goes to zero as 
n/i goes to infinity, and the frequency approaches P _1 as n/£ tends to zero. Furthermore, 
it follows from equation (|8|) that the frequencies of modes of large degree are function of 
n/£, and that the frequency spacing between modes of adjacent n for a fixed £, 5u n e, is 
approximately equal to u 2 £ P Q /£. The total number of distinct frequencies v n ni for which 
both n' and £' are of order £ scales as I 2 . These frequencies are distributed between 
and Nb- Statistically therefore, given a typical primary frequency v a < 2N B , the lowest I 



such that \v n i — {y a /2)\ < 5u should be £ ~ £ min ~ \jv a j2bv. It is worth noting that the 
frequency of the 1 = 2 f-mode (strongly excited during tidal capture) is in fact less than 
2Nb in Sun-like stars, so resonant g-modes will exist. 

The radiative damping time of long-wavelength g-modes is of the order of the photon 
diffusion time or thermal time t t h, which for a solar type star is about 10 7 years, and it 
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decreases as £~ 2 with increasing mode degree. (The mode damping time for a completely 
radiative star decreases much more rapidly with £ since the outer turning point of g-modes 
moves outward with increasing degree. However, main sequence globular cluster stars have 
a convective envelope of thickness ^ 10% of the stellar radius, and for these stars the radius 
of the outer turning point of high-degree g-modes is independent of £ at fixed frequency, 
and the damping time scales as £~ 2 ). Thus the damping rate of a mode of degree i (and 
n ~ £) is 

It ~ — • (10) 

Hh 

Parametric instability can exist only if there are modes that have both 7^ and |Au;| less 
than [see eq. (H)]. This dual condition places a lower limit to the value of 

A) K a/3i/3 2 a t which parametric instability occurs, which we now calculate. By the arguments 
of the last paragraph, the minimum degree at which the resonance condition (|7J) is met is 

/ 1 \ 1/2 

(11) 



A K a 



) 

The maximum degree is set by the requirement that mode damping be less than 
3o; Qi AoK a /3 ll g 2 /2, which gives 

" max ~ I ^ 



(12) 



Obviously £ max > £ min in order that there exist g-modes resonantly coupled to the primary 
mode. This leads to the following condition on mode amplitude for parametric instability 
to operate: 

AoKaP ^ * (13) 

For a solar-type main sequence star, typically K aftft ~ 10~ 24 erg s -1 / 2 , t th ~ 10 14 s, and 
u) a ~ 2 x 10~ 3 rad/s. Therefore, the f-mode must have energy ^ 10 36 erg in order to spawn 
a lower frequency mode via parametric instability. Solar p and f modes lie at least seven 
orders of magnitude below this threshold (see above). 

So far we have neglected stellar rotation. Rotation lifts the degeneracy of mode 
frequencies with respect to the spherical-harmonic order m: thus we have v n i m instead of u n £. 
This enlarges the number of modes available for resonant non-linear interactions, including 
parametric instability. For a star rotating rigidly with angular velocity f2 the frequency 
splitting for high degree g-modes in the rotating frame of the star, 5z/ rot = v n ^ — v n ,i-ii is 
approximately Qq/2tt£, whereas if the mean differential rotation in the radiative interior 
of the star is AQ, then 5v TO t ~ £AQ/27r (see Appendix B). Thus for a rigidly rotating 
star rotational splitting becomes important at mode degree £ such that 5u TOt ~ v/£ 2 , or 
£ w 27rz//f2 ~ 10 3 for a star like the Sun. For £ greater than this value one can find modes 
with frequency within Q / ' £ 2 of some desired value, which is much better than what we get 
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for a nonrotating star. Note that even a very small differential rotation in radiative interior 
leads to a large value of 8v rot . However, in the case of the Sun, where the internal rotation 
is known from measurements of p-mode frequency splitting, the differential rotation in the 
radiative interior appears to be close to zero. Nevertheless, in rotating stars one can find 
modes with smaller values of |Ao>| and so parametric instability persists down to smaller 
mode energies. 



2.2. Energy transfer to small-scale modes 

The astrophysical importance of the parametric instability depends upon the amplitude 
at which nonlinearities intervene to stop exponential growth. In simple mechanical systems, 
it is common for parametric instabilities to saturate because the natural frequency of 
the daughter mode is a function of amplitude. Although the same mechanism applies in 
principle to g modes, it does not affect the linear growth rate until the energy in a single 
daughter mode is comparable to that in the primary mode, as discussed in more detail in 
Appendix C. Other nonlinearities become important at much smaller amplitudes. 

As the energy of the primary mode decreases to balance the increasing energy in all of 
its parametrically unstable daughters, equations (||) and @ remain true if we replace the 
initial amplitude (A ) with its current value (A a ). Hence both the width of the unstable 
resonance and the peak growth rate decrease with A a . Initially unstable daughters become 
stable. It can be shown that parametric resonance involving a single daughter mode, where 
Pi = /3 2 , leads to periodic transfer of energy back and forth between the daughter and the 
primary. 

We now discuss the generalization of this result to the more realistic case where a 
large number of modes (/3) of angular frequency ujp « u a /2 are coupled to the mode a 
of frequency u a . At present, however, we ignore direct interactions among the daughters 
themselves, as well as parametric decay of the daughters into "granddaughters." For clarity 
we write down the equations for modes a and P below: 

3 

A a = i-u a J2 K apf3Bpexp(-iAujpt), (14) 
1 p 

Bp + jpBp = i3upK al3 pA a Bp exp(iAupt) (15) 

where 

Atop = u) a — 2up. (16) 
It is straightforward to show from these equations that in the absence of damping the 
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sum of the energy in all modes is conserved, i.e.0 

Furthermore, it can be shown that nonlinear interactions lead eventually to equipartition of 
energy among all the coupled modes. Thus the steady state energy of modes is 

(E a ) = (Ep) « Eq/N, (18) 

where N is the number of modes involved in parametric resonance with mode a. Initially, 
when all the energy is in mode a, daughter modes grow exponentially, and the characteristic 
time for the energy in mode a to halve is of order rj' 1 [eq. (El)]. We shall call this 
energy-halving timescale t n \, since it derives from nonlinear mode couplings. 

In more detail, t n \ depends logarithmically on the mean initial energy per daughter 
mode, Epfl. The latter energy has at least the thermal value kT ~ erg, whence 
t n \ < \n(Eo/NkT)/2r] « 50/r). Since the g-modes may be stochastically excited by 
macroscopic processes such as semiconvection in the core or convection in the envelope of 

I/O 

the star, it is probably more realistic to assume t n \ ~ IO/77, or approximately 50E a ^ 2 days. 

In the simplified model of eqs. (|14|)-(|T5|), complete equipartition [eq. (|I5D 1 is obtained 
only on a timescale ^> t a \, because the nonlinear interaction rate decreases in proportion 
to the primary amplitude; furthermore, phase coherence between the primary and the 
daughters tends to be lost. Because of damping and because of coupling to granddaughter 
modes, equipartition may never be achieved in a real star. 

We have verified that our model behaves as described by numerically integrating the 
differential equations (p^-fll^f), and the results are shown in fig. 1. 



We can compare the nonlinear timescale to the timescale for decay of the primary 
mode by linear radiative damping. The former is 

tnl-~- o A ~ —A • ( 19 ) 

Tj 6UJ a AoK a p 1 p 2 Va-n-QKotPxPi 

If the radiative damping time for the primary mode is l/j a , then parametric instability is 
the dominant damping process for the mode provided that 

27a 

AoKafrfa £ — . (20) 



3 Strictly speaking, the conserved energy includes terms of cubic and higher order in the 
amplitudes. In a hamiltonian formalism, the latter terms are responsible for mode coupling 
and for transfer of energy among resonant linear modes. For a weakly nonlinear system, 

however, terms beyond the quadratic can be neglected when evaluating the total energy. 
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The damping time for solar f and low order g modes is about 10 4 years. By this comparison, 
nonlinear damping would appear to be more rapid than linear damping if the energy of 
the primary is greater than 10 34 erg, which is less than the minimum energy at which 
parametric instability can exist [eq. [13])]. So the nonlinear process should always dominate 
whenever it is allowed at all. But we shall revise this conclusion below. 

We now consider the long-term evolution of the system. After a period ~ i m , the energy 
in the primary mode drops by a factor ^ 2 and is shared among iV ~ £ max (Au;/u;) ~ ^ m ax/^mm 
parametrically resonant daughter modes. Using eqs. (|TX[) - (|T2D and the typical values of 
K ai3!f3 2 , tth-, and u a cited above, we estimate N ~ 10 10 Eq^ 2 . The time scale for further 
energy decline is set by the daughter modes, which both dissipate energy by radiative 
process and also spawn lower frequency modes. The radiative damping time for g-modes of 
i ~ 10 4 is on order a week; and in fact this time is ~ r]~ l for modes of i ~ 4nax- Energy 
dissipated this way can clearly never return to global oscillations of the star. 

Additionally, energy can also be removed from daughter modes if they themselves 
parametrically destabilize modes of still lower frequency — "granddaughter modes". The 
parametric growth rate of a granddaughter mode (7) due to nonlinearity in the daughter 
mode (/3) is rf « uipKp^Ap. Energy transfer to the granddaughters will significantly retard 
the growth of the daughter modes when rf > i], or in other words Kp^Ap > K a ppA . 
Now it follows from an analysis of the 3- mode couplings that Kp 77 ~ (£ p/£ a )^ a pp if 
£/3 ^ 4* ~ 2 (Appendix A). Hence secondary parametric instabilities become important 
when A/3 ~ A /£p. Most of the daughter modes have degree £p ~ £ m ax, and after a time 
~ t n i each has a typical amplitude Ap ~ N~ 1 I 2 Aq ~ (^ m in/^max) 1 ^ 2 ^o- This is smaller than 
V4>ax if C > Cn- From eqs. flD])-© we find W4in « 25£ 3/ 4 4 2 - 

Thus in the regime E 10 40 erg, most of the daughter modes will be limited 
to energies much less than the equipartition value (Eq/N) by secondary parametric 
instabilities. We can estimate how much this process should increase the nonlinear damping 
time of the primary mode. We assume that each daughter mode (/?) continues to absorb 
energy from the primary at a rate ~ 2rjEp even after Ep has saturated. If saturation occurs 
at Ep ~ (A /£p) 2 , the primary loses energy at the rate 



~vE a , (21) 





instead of E a = —i]E a as for E a ^> 10 40 erg. The modified nonlinear damping rate implied 
by eq. fl2"T|) scales as E^J 4 instead of E]j 2 . Comparing this modified damping rate to the 
linear radiative damping of the primary mode (~ 10 4 yr), we find that nonlinear damping 
continues to dominate provided E a ^ 10 37 erg, which is only an order of magnitude above 
the threshold for parametric instability (§2.1). 

At the energies expected of tidal capture (Eq ~ io 42-45 erg), the average energy per 

5 /4 

daughter mode is so small, because their number N cx E is so large, that the growth 



-10- 



rate of the granddaughter modes can be neglected for the most part. However, secondary 
instabilities will prevent the fastest-growing daughters — those near the center of the 
instability strip (0) — from exceeding an amplitude ~ A^/lp. By contrast, in the simplified 
model of eqs. (|T^-(|l"5|) where secondary instabilities are ignored, the single fastest-growing 
mode briefly attains almost the entire energy of the system. 

Low order g-modes of a radiative star carry a substantial amount of energy as a result 
of tidal interaction. The parametric instability described above will operate on these modes 
just as well, and independently of the f-modes, causing a loss of their energy on a time scale 
similar to the f-mode. This is discussed in more detail in 84. 



3. Nonlinear energy dissipation in convective stars 



The parametric nonlinear interaction mechanism described in the last section does not 
apply to fully convective stars, which do not have stable g-modes. However, even for these 
stars nonlinear mode coupling can be an important dissipation mechanism for the f-mode 
energy. The energy in f-modes of a convective star after tidal capture is often an order 
of magnitude larger compared to a radiative star of same mass and orbital parameters, 
and thus the nonlinear parameter, A a K a ^^ is larger for convective stars. In this case even 
the normally very inefficient process of nonlinear coupling of a f-mode to higher frequency 
p-modes can transfer energy to short wavelength waves where it is dissipated on a relatively 
short time scale. The nonlinear response to the f-mode can be treated as a forced excitation 
of all p-modes of degree I = 0, 2, and 4, in which the forcing function is quadratic in the 
f-mode amplitude and oscillates at twice the f-mode frequency. We calculate this response 
by solving an inhomogenous wave equation for the sum of the p-mode amplitudes with the 
nonlinear forcing function as the inhomogenous term. The latter can be obtained either 
from the Lagrangian given in Appendix A, or by expanding mass, momentum and energy 
equations to second order in wave quantities [e.g. pzimbowski (1982)| |. The solution of 



the inhomogeneous equation is shown in figure 2 together with the f-mode eigenfunction. 
The turbulent viscous, and thermal, dissipation of wave energy is calculated using this 
wavefunction. Note the generation of significant power at short wavelengths and the 
relatively large displacement amplitude and its derivative in the outer part of the star; these 
effects enhance the dissipation of f-mode energy. The resulting nonlinear dissipation time 
for the f-mode by coupling to p-modes is 

t ffp = 3 x lO 4 ^ 1 days, (22) 

where is f-mode energy in units of 10 42 erg. 

We can estimate the damping rate in another way that is less accurate but provides 
some insight. In the present special case of the general 3- mode interaction P), we drop all 
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nonlinear terms except the square of the f-mode amplitude. The result is 

A/3 + uj 2 p A p = -3u 2 pKp aa A 2 a , (23) 

where A a is the amplitude of the f-mode which we take to be A sm(uj a t). Substituting this 
into the above equation we find for the amplitude of mode /3, 

2>ujlAlKj3 aa cos(2uj a t) 

A * = in-*® ■ (24) 

The net displacement associated with the p-mode amplitudes Ap is 

£=T,Mfi> ( 25 ) 

P 

— * 

which can be used to calculate the dampings. The function Ap£p has alternating signs for 
even and odd order modes near the surface of the star where most of the damping is taking 
place. This can lead to significant cancellation in the total displacement. Nevertheless we 
proceed by assuming that different modes can be added in quadrature. The energy in mode 
P is 



{Ep) 



9 Up(A Kp aa ) Aq 
16 



1 + 



<4 j 



(26) 



If the rate of dissipation of energy of a mode (3 driven at frequency 2u a is Tp, then the rate 
of energy dissipation of the f-mode (mode a) can be calculated from this last equation, and 
the resulting dissipation time is 

tffp-Ul^TpK*) . (27) 

The 3-mode coupling coefficients, ft// p , involving a quadrupole f and a radial p mode for 
a solar model, have a typical value 10~ 24 erg -1 / 2 , and the damping time for high order 
p-modes (l/r p ) is about 10 days. Since there are about 100 p-modes of degree 0, 2, and 4 
coupled to the quadrupole f-mode, we see from equation ( P7[ ) that for a f-mode of energy 
10 42 erg, the non-linear damping time is about 10 4 days which is within a factor of a few of 
the numerically calculated value given in eq. Q2"2]). For comparison we note that the linear 



turbulent viscous damping time of f-mode is about 10 years. 

The energy in the f-mode of a 0.5 solar mass polytropic star of index 1.5, tidally 
captured by a one solar mass star, is about 10 45 erg ps 0.1), and so its nonlinear 

dissipation p2|) j time is about 30 days. At such high energies, however, the system is very 
nonlinear: that is, four mode couplings are as important as the three mode couplings. 
Figure 2 shows that for this energy the amplitude of the p-mode response is about the same 
(or greater) as the f-mode amplitude throughout the star, and thus higher order couplings 
must be equally important. In this case it is very likely that high order p-modes are 
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resonantly excited leading to even more efficient dissipation to f-mode energy. Moreover, 
since most of the energy is dissipated in the outermost envelope of the star at a rate of 
~ 10 39 erg s _1 , which is more than can be transported by convection, it will cause the 
envelope to expand away from the star. This is a further indication of the breakdown of the 
weak nonlinear calculation presented here. 

In addition to the nonlinear effect arising from the third order fluid Lagrangian, 
considered here, there is also nonlinearity associated with the thermal dissipation process 
itself. It can be shown that both of these are of the same order, and thus a neglect of 
the latter does not change our result by more than a factor of a few, which is within the 
uncertainty of the calculation of this section. 



4. Application to capture binaries in globular clusters 



Main sequence stars of mass greater than 0.5 Mq have a sizable radiative core (see 
table 2) that supports a dense family of g-modes. When such stars are tidally captured, 
resulting in a highly eccentric orbit with separation between the stars at periastron of 
only a few stellar radii initially, their quadrupole f and low order g modes are excited 
to large amplitudes. The linear radiative damping time for these modes is of order 10 4 
years. However these modes are resonantly coupled to high degree g-modes by parametric 
instability which provides a very efficient way of dissipating their energy (see §2). We 
now estimate the time for the energy of the low order modes to decay due to nonlinear 
interaction using the results of §2. 

The main ingradient required for nonlinear dissipation calculation is the 3-mode 
coupling coefficients, k,,^^, which are computed using equation (|A17|) . The numerical 



calculation of Kafiifc is carried out for a simplified polytropic stellar model, which is 
widely used in tidal capture and evolution calculations, as well as a standard solar model. 
The main sequence globular cluster stars of mass greater than about 0.5 Mq have a 
radiative core (see table 2) and can be crudely modelled as a polytrope of index 3 [cf. 
Pinsonneault, Deliyannis, fc Demarque ( 1992 )| , [Lai, Rasio, fc Shapiro (1994)| , Peliyanms| 



1995)|| . The coupling coefficients involving low degree g-modes (lp < 40) are calculated 
using eigenfunctions obtained by solving the wave equation, whereas for high degree 
g-modes we use the WKB solution which is an excellent approximation. We find that 
Ka/3/3 for the Sun and 1.0 Mq polytropic star of index 3 agree to within a factor of two 
f° r £p ^ 100 (the polytropic model yields a higher value). However, coupling coefficients 
increase with £p for polytropic models but not in the case of the Sun or a realistic stellar 
model. This is because g-modes are trapped in the radiative interior of low mass main 
sequence stars, whereas in a polytropic model, where the Brunt- Vaisala frequency increases 
monotonically with distance from the center, g-modes propagate closer and closer to the 
surface with increasing mode degree. One should therefore use a realistic stellar model 
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for the purpose of nonlinear mode interaction calculation. Since K a p x fo ~ 
depends only weakly on stellar mass we use a solar model to compute coupling coefficients 
that are used for lower mass globular cluster stars; the true coupling coefficients are larger 
by about 20%. Table 1 shows the value of K a pp for a selected number of modes. Note that 
a typical value of coupling coefficient is 10 -24 erg -1 / 2 which is in good agreement with the 
rough estimate of G~ 1 I 2 B^J 2 M~ X . 

The dependence of n a pp on mp is entirely contained in Wigner 3-J symbol [see eq. 
( [A10|) l, which for large tp scales as cos^mpfelp). Thus a good fraction of modes of all 
possible m, for a given £, are equally well coupled to the f and low order g modes. Finally, 
we note that the coupling coefficient decreases with increasing An = \n^ — rip 2 \ and is 
neglegible if An ^10. 

The energy input rate per orbital period, for several different modes, are shown in 
figure 3 as a function of distance between the stars at periastron. This was calculated 
using a scheme similar to |Press fc Tcukolsky (1977] . In general most of the energy is in 



some g-mode which has a period close to the periastron passage time. Substituting mode 
amplitudes and coupling coefficients into equation (|P9|) we obtain nonlinear damping times 
which are displayed in figure 4 for several modes as a function of periastron separation. 
Note that for mode energy of 10 42 erg (corresponding to dimensionless surface displacement 
amplitude 5R/R* ~ 2xl0 -3 ; see Table 1) the nonlinear damping time is 50 days. During 
this time the mode energy decreases by a large factor, which is shared equally by about 
N ~ 10 10 high degree g- modes that have been parametrically excited via nonlinear mode 
coupling. It follows from equation ( |19|) that the nonlinear damping time scales as E~ x ' 2 , 
and thus for mode energy of 10 44 erg the damping time is only 5 days. For comparison we 
note that the radiative damping time for low order quadrupole modes are of order 10 4 years. 
Thus nonlinear dissipation dominates even when periastron distance is 4R* (see fig. 3 & 4). 
In fact the nonlinear dissipation rate is larger than the linear radiative damping so long as 
mode energy is greater than about 10 37 erg. It should be emphasized that energy removed 
from a mode by nonlinear interactions is distributed over a large number of g-modes. The 
energy is eventually thermalized on the dissipation time scale of high degree g-modes, 
which for modes of £ ~ £ max is of the same order as t n \. How the star adjusts to the energy 
deposited in its core on a short time scale is outside of the scope of this work. However, we 
note that even 10 44 erg of energy deposited in the radiative core is a small fraction of its 
thermal energy content, and thus will have a minor effect on the structure of the star. 

Main sequence globular cluster stars of mass less than about 0.5 M Q are almost fully 
convective, and therefore the coupling of f-mode to g-modes is very weak. The dominant 
nonlinear coupling of the f-mode in this CclSC, clS described in §3, is with p-modes. Physically 
this nonresonant mode coupling results in power on smaller length scale, and also the value 
of the gradient of the wave function generated can be considerably larger than the original 
f-mode wave function in the outer part of the star (see fig. 2). This leads to more efficient 
dissipation of wave energy. As discussed in §3 the f-mode dissipation time due to this 
process is about 3xl0 4 Ej 2 days. 



-14- 



5. Summary 

Tidally excited modes of a main sequence star in a highly eccentric orbit, with 
periastron distance of 4 stellar radii or less, have energy about 10 42 erg or greater, or 
fractional change in stellar radius of 0.2% or more. The linear radiative and turbulent 
dissipation times for these modes are about 10 4 years. We have calculated the time needed 
for energy of tidally excited modes to decrease due to nonlinear interaction with other 
modes in the star for two different cases that are described below. 

One of the cases considered is that of a star with a radiative core - for instance 
main sequence globular cluster stars of mass greater than about 0.5 Mq. The most 
important mode coupling in these radiative stars involves high degree g-modes which are 
parametrically excited by the low order quadrupole modes. The time scale for energy 
transfer for this process is 2/ ' {y a K a ppEH 2 ), where K a pp ~ G~ X I 2 M~ X R X J 2 ~ 10 -24 erg -1 / 2 
is 3-mode coupling coefficient. Thus for E a m 10 42 erg, or 5R/R* ~ 2xl0 -3 , the nonlinear 
damping time is 50 days. This is smaller by a factor of 10 6 compared to the nonlinear 
damping time estimated by Kochanek (1992), who had not considered parametric coupling. 
The energy removed from the primary mode initially resides in roughly 10 10 high degree 
g-modes and is later dissipated in the radiative core of the star on a time scale of about a 
year. Since this energy is small compared to the thermal energy content of the core, it is 
unlikely to have much effect on the structure of the star in one orbital passage. 

The other case we have considered is that of a fully convective star, such as a main 
sequence star of mass less than about 0.5 Mq. Most of the tidal energy in these stars, 
which do not have stable g-modes, resides in quadrupole f-modes. The f-modes couple 
to high order p-modes of degrees £ =0, 2, and 4 generating short wavelength waves that 
significantly enhance the energy dissipation rate. The energy of short wavelength p-modes 
is dissipated near the stellar surface by turbulent and radiative viscosity. The nonlinear 
energy dissipation time of f-modes is estimated to be 3 x 10 4 (10 42 erg /Ef) days or 2 x 10 _1 
(5R/R*)~ 2 days, which is shorter than the linear dissipation time when the mode energy 
is greater than about 10 40 erg. And it is also smaller than the nonlinear damping time 
estimate of Kochanek (1992) by a factor of about 10 3 . For mode energy of 10 45 erg, that 
which is expected in tidal captures, the nonlinear damping time is about 30 days. However, 
for such a large energy, 4- mode 5-mode, and higher-order mode coupling are also important; 
that is, the weak nonlinear approximation is not valid. 
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A. Adiabatic g-mode couplings in spherical geometry 

The Lagrangian for inviscid, compressible fluid with internal energy E(s,p) is 

C = p(j-E{8,p)-U{x)y (Al) 

where v is fluid velocity, s, p, and U are entropy, density, and gravitational potential. 
The three mode coupling is described by the fluid Lagrangian expanded to third order in 
displacement and is given below 

c 3 = -^e'% ; ,-%rv-e-^(v-i5> 
+ ^p^(v • a 3 - ~p ; *e(v • a 2 - v • a 

- \p;i;y,k€^e ~ \u.^, h p f£ J ~ \pC^SU^. (A2) 

Here p(r) is the equilibrium pressure profile, £ is the fluid displacement from its equilibrium 
position, semicolon denotes covariant derivative (see Misner, Thorne, and Wheeler 1973 for 
the definition of covariant derivative), and 5p is the Eulerian pressure perturbation to first 
order in £: 

5p = ■ Vp - — V ■ £• (A3) 

Let £ — > £ + S, where H is a large-scale fundamental or g-mode and £ becomes the sum 
of all small-scale modes. Then the dominant (i.e., leading order in i of small-scale modes) 
terms are 

+ ^^^Uv • a 2 v • s - ■ if-p-.cv ■ £V • s 

- /':-:, Qaev • s + a^v ■ a) - ^(p*** + m 

Capital letters are used to distinguish quantities pertaining to the large-scale mode: e.g. 
5P and 5U are the first-order eulerian pressure and potential perturbations associated 
with H. (However U unadorned with 5 is the equilibrium potential.) The 3-mode coupling 
coefficient is the integral of £3 over the star, when the mode eigenfunctions are normalized 
to unit energy. 

Integration by parts has been used in deriving formulae (|A2|) and (|A4| ) to remove as 
many derivatives as possible from the small-scale displacements £, except where these occur 
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as divergences V • £. In general ~ £/A, where A <C .R* is lengthscale associated with 
the displacement. Because the small-scale g-modes are nearly incompressible and isobaric, 
however, V ■ £ ~ £/H and 5p/p ~ C.X/H 2 , where if ~ i2* is the pressure or density scale 
height. 

The next step is to integrate £ 3 over solid angle (d 2 Q = sin 9d9d<p) . As usual, the 
covariant components of the angular displacements can be written as partial derivatives of 
a scalar potential [here and henceforth, greek indices range over the angular coordinates 
{9, 0}, whereas roman ones range over {r, 9, 0}], 

f a = d^, ~ a = d a y, (A5) 
and the displacements can be decomposed into angular harmonics: 

C = E&(*>r)n(M), ^EW.W.A (A6) 

A A 

E r = £~^,r)y A/ (M), * = J2* Af (t,r)Y Af (9,<f>). (A7) 

A/ A/ 

We use the single symbol A for the degree and order (£, m) of a spherical harmonic i.e., 
Y A = Yg m (9, (ft), and define A 2 = £(£ + 1). The projection of a function f(t, r, 9, (p) onto Y im 
is written f A and is understood to be a function of (t, r) (except that Y A denotes Yg m (9, <j>)): 

f A = J f(t,r,9,<P)Yl(9^)d 2 n. (A8) 

Enclosing a term in angle brackets means to integrate it over solid angle; thus the above 
could be written f A = (Y A f). The angular integral of three spherical harmonics is related 
to the Wigner 3j symbol (see Edmonds, 1960), 



/ *1 *2 *3 
Y 777.1 m 2 m 3 

by 

(AiAsAg) = (lA^A^Aa) 

~ (2li + l)(2l 2 + l)(2l 3 + P n 1/2 

47T 



(A9) 



«1 «2 «3 «1 *2 «3 

mi m 2 m 3 / I 



.(A10) 



To calculate (£ 3 ) from eq. (|A4j) , one needs the following identities: 

,2 , A2 _ A2 



Y A:i g aP (V a Y Al )(VfsY A2 )) = ^i±^^(AiA 2 A 3 ), (All) 



V^(v,v,r A3 )(v a r Al )(v^ A2 )) = ( A ') 2 ( A i A i) 2 (AiA2As)- (A12) 
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Here V is the usual covariant derivative on the sphere; that is, V 7 g a( g = 0, where g a/ 3 is the 
2D metric 

g a pdx a dx p = dd 2 + sin 2 ddcj) 2 , (A13) 

and g a @ is its inverse. 

To establish Eq. ([All ), for example, note that 

r^(v a y Al )(v^y A2 ) = \ [v 2 (r Al r A2 ) - r Al v 2 y Aa - (v 2 r Al )r A2 ] , (am) 

and that 

V 2 F A = g a/3 V a VpY A = A 2 Y A , (A15) 

and use integration by parts. The 2D covariant derivative is not simply the restriction of 
the 3D covariant derivative to the sphere, but is related to it by formulae such as 

f;a = V a / = d a f, 

r t~7 v — 1 i r — 1 

V -a = VaV - T V a , O a V - T V a , 

v a; p = VpV a - rg a pv r , (A16) 

where / and v 1 are scalar and vector fields in 3D. To avoid confusion, we raise and lower 
indices using only: for example, v a = g a iV l (not g a p^)- 

Using these rules we integrate the third order Lagrangian over angle which leads to the 
following expression 



/ C 3 d 2 n = J2 (AiA 2 A/> { 

J Ai,A 2 ,A/ I 



2\2 



(A/) 2 - (Ai - A 2 2 ) 
4r 4 

V + A^-A 2 r 

2r 2 



* A/ -0 A i<*P A 2 - ^ ( 5P + P SU ) Af ^A!^A 2 



+ 



if 1 dp -> - 

rE r A ip Al SpA 2 - - rd r 8P + prd r 6U + --V • S 
1 2 \ r dr 



+ 



d ( 1 dp\ _ r 
dr \rdrl 



+ 



A/ + Al 2 _ A2 2 
2r 2 
-l. 



vS r - ^) Af ip Al d r 5p A2 + rd r (^) Af 8p Ai a 2 - rd r (^5P) Af ij Al C A2 



- P rd r (r^5U) A ^ Ai e A2 - [ SP + j ^ ^ V • £ Aa - ± [ ^ ) ^A^ 



d / 1 dp N 
(ir \r dr 



+ 



a r H A/ e Al a r 5 PA2 - (d r 5P + • h + |§h^ e Al v • 4 
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~ Up + 7(7 - 2) P V • 3 + ^E r ) V • f Al v • d 2 



(A17) 



Every term in this expression is formally of order (£ r /H) 2 p, where £ r = £ A or £ A , since 
t/> Ai ~ ^Xi/^- The 3-mode coupling coefficent, « a( g 7 , is obtained by integrating the above 
expression over the radius of the star i.e. 



K af 3y = J drr 2 J d 2 QC 3 , (A18) 
Note that all the eigenfuctions in equation (A17) are normalized to unit energy. 

When all three interacting modes have wavelengths (Ai,A 2 ,A 3 ) <C then the 

first term on the righthand side of (|A2| ) dominates the rest by a factor ~ H/ max(Ai, A2, A3). 
This justifies our statement in the text that the coefficient K/3 77 coupling daughter to 
granddaughter modes exceeds the coupling K a pp of the primary mode to its daughters by a 
factor ~ I p. Physically this can be understood, or at least remembered, by noting that the 
parametric growth rate scales as the shear ~ ou^/X in the parent mode. 
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B. Rotational splitting of eigenfrequencies 



The rotational splitting of mode frequencies as seen in the mean rotational frame of 
the star is given by (e.g. [Stix (199TJ ) 



A^ m = v nlm ~ v nlQ — —mC n £ri /2lT, (Bl) 

where Qq is the mean angular speed of the star, 

r I d 3 x p [ti + 2Uh - (Ml/Op) - q) 2 + (f + £ - 2)ti}} 

and 

Afi(r) = Q(r) - Q - (B3) 

For high order g-modes £ r ~ i££,h, and thus / d 3 x p£ r £h ~ 0. Substituting these into 
equation (|B2| ) we find that for a rigidly rotating star 

Cm « ^, (B4) 

and so the frequency difference between modes of m = i and — i is Qq/27t£. For a 
differentially rotating star the splitting coefficient is 



Cnt w 2^ ~ T^' (B5) 



where Af2 is the mean value of differential rotation in the radiative interior of star where 
high order g-modes are trapped. This leads to Ai/ n £ m « — mAf2/27r if £ ^> 1. However, if 
this formula predicts Az/ n ^ m > N B , the Brunt-Vaisala frequency, then individual g-modes 
are trapped within regions smaller than the entire radiative zone, and the spectrum of mode 
frequencies becomes continuous. 
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C. Frequency shifts and saturation 



Mechanical oscillations and waves are generally anharmonic: their natural frequency 
depends upon amplitude. Unless other processes intervene earlier, a parametric instability 
will saturate at an amplitude such that the resonance condition (^j) is no longer satisfied 
because of the anharmonic frequency shift. 

For example, a generic one- dimensional oscillator with Hamiltonian 

H = ^ + ^fi_ + " + » q > +olitM \ (C1) 



has free oscillations (cf. [Landau fc Lifshitz 1976| ) 



q(t) = a cos[uj(a)t] 
10(a) = io 



12 J \a 



(C2) 



Here ao, a characteristic amplitude at which anharmonicity becomes important, is chosen 
so that the larger of p and k? is of order unity. Suppose this oscillator is coupled to another 
in such a way that very small oscillations are parametrically unstable with growth rate 
7] = eujQ. The range of frequencies over which parametric instability occurs is Au « rj. If 
(for simplicity) parametric resonance is exact at infinitesimal amplitude, then growth stops 
when \uj(a) — luq\ « rj, hence a ~ e 1 ^ 2 aQ. 

Waves in fluids are normally strongly nonlinear when ka ~ 1, where a = max(|£|) 
is the peak displacement and k = 27r/A is the wavenumber. That is, we may normally 
take ao = X/2tt. This is because the eulerian dynamical equations contain the advective 
derivative d t + v ■ V, and for small perturbations of a static equilibrium, the ratio of the 
second (nonlinear) part of this operator to the first is 0(kSv/u) = 0(ka). A single g 
mode, however, is an exception to this general rule if the wavelength is small compared 
to a pressure or density scale height, H. Such a wave is nearly incompressible, so that if 
Sv oc exp(ik-x) in a local plane-wave approximation, then V-ot/is ~ 5v/H instead of ~ kSv. 
It follows that the operator v ■ V = 0(Sv/H) when applied to any fluid perturbation whose 
WKBJ wavenumber is parallel to that of Sv. One can check that other nonlinear terms 
vanish to leading order in (kH)' 1 . For example, the dynamical equation for the vorticity 
V x v contains the source term p~ 2 (Vp x Vp); but for a single wave, the gradients of the 
density and pressure perturbations are both approximately parallel to the wavenumber, so 
that VSp x VSp vanishes to leading order in (kH)^ 1 . 

We note in passing that nonlinear interactions involving two or more nonparallel 
wavenumbers need not be similarly suppressed. In parametric instability, for example, the 
primary and daughter modes have nonparallel wavenumbers. Anharmonicity is a special 
nonlinearity because it involves self-interaction of a single mode. 



Consequently, the anharmonic frequency shift of a WKB g mode is 



co(a) — ujq — O I ujq 




) 



(C3) 



Notice that this does not depend on mode wavelength. Equating this shift to the width 
Auj ~ eujq of the parametric resonance caused by a fundamental mode of dimensionless 
amplitude e = dR^/R*, one predicts saturation when the small-scale mode has amplitude 
a ~ e 1 / 2 H. In the radiative core, however, H ~ R, so the small-scale mode would have 
more energy than the large-scale one by a factor of order (a/5R*) 2 ~ e _1 . 

Clearly, other nonlinear processes must intervene to halt the growth of the small-scale 
mode before frequency shifts become important. In the situation contemplated here, we 
foresee only two such processes: one is the decrease in the parametric growth rate as the 
energy in the primary mode declines; and the other is the destabilization of granddaughter 
modes. Both processes have been discussed in §2.2. 
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Table 1 

Mode coupling coefficients for quadrupole f- and g-modes of a solar model with high degree 
g-modes of half the frequency. First two columns are the number of radial nodes and the 
frequency of the primary mode (a). Third column is coupling coefficent for parametric 
decay into identical daughter modes Last column gives the fractional surface 
displacement amplitude of the primary mode when its energy is 10 42 erg. 
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Table 2 

The size of convection zone in main sequence stars of metalicity (z) equal to 0.0001 and age 
13.0 Gyr. The first column is the mass of the star in the unit of solar mass, the second 
column is the effective temperature, the third column is the thickness of the convection 
zone divided by the radius of the star, and the last column is the fractional mass in the 
convection zone. These numbers, obtained using the Yale stellar evolution code, were 
kindly provided by Constantine Deliyannis. 
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Figure Captions 

FIG. 1. — Parametric nonlinear interaction involving 10 3 modes. The amplitude of the 
primary mode is shown as a function of time. 

FIG 2. — The upper panel shows the radial displacement eigenfunction of the f-mode of 
energy = 10 45 erg (solid line), and the superposition of radial p-mode displacement-functions 
generated as a result of nonlinear couplings to the f-mode (dash-dot curve). The lower 
panel shows the radial derivative of wave functions shown in the upper panel. 

FIG. 3. — Energy input per orbit (in erg) into quadrupole f- and a few low order g-modes 
of m = £ = 2, as a function of periastron distance. The mass of the star is 0.7 Mq, and 
the orbital eccentricity is 1.0. The star has been modelled as a polytrope of index 3 with 
the ratio of specific heat taken to be 5/3. From top to bottom at left (d per i/-R* = 2.5) the 
modes are g 2 , /, g 4 , 9e, 9s, and g 10 . 

FIG 4. — The log of nonlinear damping time (in days) for a few low order quadrupole 
modes of a 0.7 Mq radiative star as a function of periastron distance. The linear radiative 
damping time for all of these modes are about 10 4 years. Curves correspond to same modes 
as in Fig. 3: that is, from top to bottom at d peT i/ R* = 5, /, g 2 , #4, g%, gs, and g w . 



